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Abstract 

We study suspensions of solid particles in a viscous incompressible fluid in the presence of 
highly oscillatory velocity-dependent surface forces. The flow at a small Reynolds number 
is modeled by the Stokes equations coupled with the motion of rigid particles arranged in 
a periodic array. The objective is to perform homogenization for the given suspension and 
obtain an equivalent description of a homogeneous (effective) medium, the macroscopic 
effect of the interfacial forces and the effective viscosity are determined using the analysis 
on a periodicity cell. In particular, the solutions u^ to a family of problems corresponding 
to the size of microstructure e and describing suspensions of rigid particles with random 
surface forces imposed on the interface, converge H 1 - weakly as e — > a.s. to a solution of 
the so-called homogenized problem with constant coefficients. It is also shown that there 
is a corrector to a homogenized solution that yields a strong H 1 - convergence. The main 
-^i . technical construct is built upon the T- convergence theory. 

Key words: effective viscosity, velocity-dependent forces, homogenization, Stokes equa- 
t^j- ' tion 
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1 Introduction 

Flows of incompressible fluids that carry rigid particles also known as particulate flows 
are essential parts of many engineering and environmental processes (i.e. particle sedimen- 
tation, fluidization and conveying) and are commonly encountered in many applications and 
fundamental fluid mechanics. The complexity of mechanisms that govern fluid-particle and 
particle-particle interactions makes the numerical simulation of these flows be one of the most 
challenging problems in computational fluid dynamics. Many analytical, numerical, and ex- 
perimental studies have been performed during the past decades, however still much more 
research is needed for the fundamental understanding of these complex heterogeneous media. 

For instance, the Theological behavior of particulate flows have been extensively studied 
for over a century. Up to date the most investigated concentration regime of particles is the 
dilute one in which the hydrodynamic interactions between inclusions are normally ignored 
as the interparticle distance exceeds the range of the flows that appear due to the particle 
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motion. Hence, it is possible in this case to isolate the effect of particle-fluid interactions on 
the effective behavior of suspensions, as was done in the pioneering work by Einstein [15] where 
the asymptotic study of the effective viscosity n* of a suspension of rigid particles in a small 
concentration regime was carried out. A simple approximation to the viscosity \x* in volume 
fraction <fi of rigid neutrally buoyant spheres in the suspension had been formally derived there 
(see also [16]): 
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M *~/m + -0) as 0^0, (1.1) 

where ^u is the viscosity of the ambient fluid. The subsequent extension to the ellipsoidal 
particles was done in [22]. Despite the seemingly simple linear relation (1.1) between [i* and 
(/), the rigorous justification of the Einstein's formula (1.1) has been carried out over a century 
later, in [20] by constructing explicit upper and lower bounds on the effective viscosity fi*. 

A pairwise particle interaction in the dilute regime has been first taken into account in 
[5] to formally derive a O(0 2 )-correction to the Einstein's formula. Such an 0((p 2 ) -correction 
in the case of a periodic suspension of spherical fluid drops of viscosity rj — >■ oo in a fluid 
of viscosity \x was rigorously proven, also very recently, in [1] using techniques different from 
[20] based on the layer-potential approach previously developed by the authors. For a more 
detailed description of asymptotic studies for the effective viscosity of dilute suspensions we 
refer a reader to [1] and references therein. 

Another "extreme" regime of particle concentration is the dense packing of particles in 
a suspension when the typical interparticle distances are much smaller than their sizes. In 
such a regime the effective rheological properties of suspensions exhibit singular behavior as 
the characteristic distance between particles tends to zero (or equivalently, as 4> tends to the 
maximal packing volume) . Such densely-packed particulate flows were also extensively studied 
both numerically and analytically (see e.g. [17, 19, 29, 32] and references therein) over the past 
decades. For example, in [17] a local formal asymptotic analysis based on a pair of closely 
spaced particles showed that the effective viscosity, described by the viscous dissipation rate, 
exhibits a blow up of order 0((5 _1 ), where 5 is the distance between the neighbors, whereas 
numerical study of [32] revealed that in some cases the blow up might be much weaker, e.g. 
of O ( | In <5 1 ) . Such a discrepancy comes from the fact that in the high packing regime the 
dynamics of particulate flows is driven by the long-range interactions between particles and 
local asymptotics is not sufficient here, therefore analysis of this concentration regime is quite 
challenging. Luckily, the development of discrete network approximation techniques of [7-9] 
allowed to settle the disagreement between the results of [17] and [32]. Namely, it was shown 
that there are multiple blow ups of the effective viscosity as 5 — > demonstrating examples of 
the cases when the stronger blow up degenerates so the weaker blow up becomes the leading 
one (see, in particular, [9, 10]). 

The third concentration regime when particles are of the same order as their sizes is called 
the finite or moderate concentration regime. Unlike in the regimes mentioned above, where 
the extreme properties in some sense facilitated the corresponding analyses (i.e. negligible 
interactions between particles in the dilute limit and the presence of strong lubrication forces 
between closely spaced particles that contributed to the blow up of the effective viscosity in 
highly concentrated suspensions), the case of finite concentrations is much harder to analyze. 
Such a regime was treated in [30] where a periodic array of spherical particles was considered. 
Under the assumption that all inclusions follow the shear motion of the fluid (formula (5) in 
[30]) it was shown that /i* = 0(e _1 ), where e is the distance between periodically distributed 
particles. This assumption is analogous to the well-known Cauchy-Born hypothesis in solid 



state physics, which is known to be not always true [18]. Also, a finite concentration regime 
was considered in [26] where an asymptotic expansion of the effective viscosity was constructed 
assuming a periodic distribution of particles of arbitrary shape. The formal two-scale homog- 
enization was carried out there under the assumption that the number of particles increases 
to infinity as their total volume remains constant. We distinguish this work [26] as being the 
closest in its goals to the current paper as stated below. 

We also point out that the definition of the effective properties of suspensions that was 
employed in [26, 30] is different from the aforementioned treatments of the dilute limit and 
dense packing regimes, where the effective properties were determined from the equality of the 
viscous dissipation rates in the suspension and the effective fluid. The latter definition of the 
effective rheological behavior is directly related to viscometric measurements that necessarily 
include boundary conditions and applied forces. 

The current paper focuses on particulate flows consisting of a viscous incompressible fluid 
that carries rigid neutrally buoyant particles whose size is comparable to the typical interpar- 
ticle distances. 

When the surface forces are negligible, the effective behavior of suspensions is well un- 
derstood. In [4, 11, 21, 24, 26] it was shown that suspensions behave like viscous fluids with a 
modified viscosity. This is also the case of emulsions studied formally in [24] and in the frame- 
work of two-scale convergence on [27]. For these systems surface tension effects are present in 
the effective viscosity. 

By contrast, other suspensions exhibit a different non-newtonian constitutive behavior. 
In colloidal suspensions the interparticle forces, including van der Waals, electrostatic, steric, 
and depletion forces, have an important role in the colloidal stability and in the suspension's 
rheological behavior. In electrorheological fluids, surface electric forces, change the viscous 
behavior in a viscoplastic or Bingham rigid-plastic behavior. In other materials electrokinetic 
phenomena and the interaction with ionized particles play an important role, as in the trans- 
port through natural clays, in the electrophysiology of cartilages and bones, in semiconductor 
transport or in membrane or bulk ion-exchangers. 

In all these examples, an important problem is the better understanding of the influence 
of the highly oscillating surface forces in the constitutive behavior of the suspension. The re- 
versible constitutive change in electrorheological fluids, the swelling in clays and the lubricating 
properties of connective tissues cannot be explained without taking into account interfacial 
phenomena and field interactions. 

The main goal of this paper is to derive and justify homogenization-type results of suspen- 
sions of rigid particles in the presence of highly oscillating random surface forces, dependent 
on the velocity field, and the microstructure in the limit of e — V 0, where e is the size of 
that microstructure. Here, we consider only noncolloidal suspensions in which hydrodynamic 
interactions are much stronger than the Brownian motions, hence, the latter are neglected. 

The relatively young homogenization theory specifically designed for analysis of highly 
heterogeneous and microstructured materials (see [3,6,23,31]) has been an active area of 
research in recent decades. In this theory the effective material properties of periodic materials 
are defined based on analysis on a periodicity cell, and then these properties depend on the 
mechanics of constituents and geometry of the periodic array of particles but independent of 
the external boundary conditions and applied forces. They are normally determined in the 
limit as the size of the microstructure e — > 0, or equivalently, as the number of particles goes to 
infinity. Such an analysis is carried out in the present paper. In particular, the solutions u^ to 



a family of problems corresponding to the size of microstructure e and describing suspensions 
of rigid particles with random surface forces imposed on the interface, converge H l - weakly 
as e — > a.s. to a solution of the so-called homogenized problem with constant coefficients. 
Convergence of the corresponding functionals describing the respecting problems is also shown. 
We also demonstrate that there is a corrector to the homogenized solution which yields a strong 
H 1 - convergence. 

Despite the fact that assumption about periodicity of a particle array is too restrictive for 
suspensions of moving particles it has been extensively used in analysis of particulate flows (see 
e.g. [17, 19, 26, 30]). Since we attempt to develop the homogenization theory for suspensions 
of rigid particles the periodicity assumption is imperative in our construct. Also, we point 
out that in the particular case when the highly oscillatory random forces on the boundaries of 
particles are absent we recover the formal asymptotics of [26]. The mathematical tools, that 
the main justification techniques are built upon, are based on Y— convergence theory [13]. 

The paper is organized as follows. In Section 2 we present fluid-particle problem setting 
including a description of interface forces and the statement of existence and uniqueness of 
the solutions etc. The homogenization problem is formulated in Section 3 and main results 
including cell problems, effective viscosity derivation, convergence and corrector results are 
given in Section 4. Conclusions are presented in Section 5. Proofs of auxiliary facts including 
existence and uniqueness of solutions are given in Appendices. 

Acknowledgments. Y. Gorb and F. Maris were supported by the NSF grant DMS-1016531. 

2 Problem Formulation for the Particulate Flow 

Let I? be a bounded domain in W* with Lipschitz boundary. We assume that D is split into 
two parts, a part Df in which we consider a viscous, incompressible fluid of viscosity /i, and a 
part D r CC D that consists of a finite union of disjoint rigid particles, D r = |J;=i D l r . Each 
rigid particle D\. is assumed to be an open connected set with Lipschitz boundary, compactly 
included in D. In the following, by n we denote the unit normal on the boundary of Df, 
directed outside the fluid region. 

The following notations are also to be used: 

d 
a ■ b = V^ dibi, for all a,b Gl" 1 
i=i 

d 
A:B= ^T AjBy, for all A, B € R dxd 

(a x b)ij = aibj - ajbi, for all a, b € R d , so a x b G R dxd . 

In the fluid region Df of the domain D, we introduce the velocity field u of the fluid and 
its pressure p. The strain rate tensor is denoted by e(u), defined componentwise by 

M«) = J, (a£ + a£J , for0<i,i<d > (2.1) 

and the stress tensor a defined by 

a = cr(u,p) = -pI + 2/j,e(u). (2.2) 



In Df, the motion of the fluid is described by the stationary Stokes equation for an incom 
pressible fluid 

- div a = f in Df , 
divit = in Df, 



(2.3) 



where / represents the body forces. On each rigid particle D r , the velocity field u satisfies a 
rigid body motion 

(2.4) 



e(u) = in dL, for each < I < N. 



We assume that for each < I < N there exists the superficial force fi = fi(x,u) acting 
on the boundary of D\. that depends on the position x as well as on the velocity field u, and 
we impose the balance of forces on every D l r by 

- / ands + / fdx + / f l (x,u)ds = 0, (2.5) 

JdD\, Jd\, JdD l r 

as well as the balance of torques 

— / an x (x — xi)ds + J f x (x — x\)dx + / fi(x, u) x (x — xi)ds = 0, (2.6) 

JdD\. JD\. JdD l r 

where x\ is the center of mass of D l r . Equations (2.3) - (2.6) are supplied with boundary 
conditions on dD that for simplicity are assumed to be no-slip conditions: 

u = on dD. (2.7) 

The boundary conditions on D l r , for each < I < N, also are chosen no-slip, thus the continuity 
of the velocity field across 8D r is imposed: 

H = on dD l r , for each < I < N, (2.8) 

where by [w] = Uf — u r we understand the jump of the velocity field u, with Uf being 
the velocity of the fluid on the boundary of D r and u r the velocity of the rigid part on the 
boundary of D r . 

Collecting equations (2.3) - (2.8), we derive the system that we intend to study: 



diver 

divii 

e(u) 

M 

crnds 



f in D f , 

in D f , 
in D r , 
on 8D r 



I 

JdD\ 
I an x (x — X[)ds 

JdDl 



= / fdx + / fi(x,u)ds, for < / < N, 

JD\. JdD l r 

= f x (x - x t )dx 
Jd\, 

+ I fi(x,u) x (x — x{)ds, for < I < N, 
JdDl 



(2.9) 



u 



IdD 

on dD 



For the system (2.9), according to the no-slip boundary conditions imposed on dD and on 
dD r we look for the solution u £ HQ(D) d satisfying divu = in Df and e(u) = in D r . Since 
the velocity field u inside each rigid particle satisfies both conditions div it = and e(u) = 0, 



the solution u is a divergence free vector field from Hq(D) . To that end, we introduce the 
following subspaces of H l (D) d 

V = {v g Hl{D) d | divi? = 0}, (2.10) 

and 

V r = {v € V | e{v) = in D r }. (2.11) 

Therefore, we look for the solution u € V r and p € L 2 (Df) satisfying the stationary Stokes 
equation (2.3) in Df in the weak sense and the equations (2.5) and (2.6) representing the 
balance of forces and torques. 

We assume that / G L 2 {D) d and for every < / < N there exists a function gi(s,z), 
gi : dD\. x R d -> R such that for all s € <9Dj. and all z £R d 

f l (s,z) = -V z g l (s 1 z), (2.12) 

where by V z we denote the gradient with respect to z, and moreover g\ : dD l r xR^E 
satisfies the following properties: 

(Gl) for every fixed z € R the function s i— > gi(s,z) is measurable with respect to the 
surface measure on dD l r , 

(G2) for every fixed s G dD l r the function z h-> ^(s, z) is convex and Frechet differentiable 
at every point z G R , and the Frechet differential is continuous on R , 

(G3) s ^ gi (s, 0) belongs to L l {dD l r ), 

2 . 

(G4) there exists 7 € (0, 1] and a function ai € L 2 -~t (dD r ), such that 

gi(s,zi) - gi(s,z 2 ) < Cai(s)\zi - z 2 | 7 , (2.13) 

for almost every s G dD l r and all Z\,z 2 € R d . 
Remark 2.1. The function g\ : dD l r x R d — > R, using (2.13) satisfies 

\9i(s,z)\ < \ 9l (s,z)-gi(s,0)\ + \g l {s,0)\ < Ccn(s)\z\~< + \ 9l {s,0)\ 
which, after using Holder's inequality becomes 

\ gi (s,z)\ < \ 9l (s, 0)1 + C^^ ai (s)^ + C^\z\ 2 

so we have the following estimate 

\gi(s,z)\<k(s) + C\z\ 2 , 
for almost every s € dD r and all z G R , where bi € L l {dD r ). 

Remark 2.2. The function gi : dD l r x R d — > R is a Caratheodory function, thus, for every 
vector field u € H l (D) d , the function, x t— y gi(x,u(x)) is well defined on dD\. and measurable 
with respect to the surface measure. Also, using Remark 2.1 the function x 1— >■ gi(x,u(x)) is 
inL l (dD l r ). 



Remark 2.3. From the convexity of z ^ gi(s,z) we have the following subdijferential type 
inequality: 

gi(s, z + w)> gi(s,z) + w ■ V z gi(s, z) 

for almost every s G dD l r and all z,w G M d , which implies that 

\w ■ V z gi{s, z)\ < 2ai{s) + C{\z\ 2 + \w\ 2 ) 

so taking all w such that \w\ < \z\ + 1 we obtain that 

\V z gi{s,z)\ <ci(s)+C\z\, 

for almost every s G dD l r and all z G M d , where c\ G L 2 (dD l r ). Also, in the same way 
as in Remark 2.2, for every vector field u G H 1 (D) d , the function x i-> V x gi(x,u(x)) is 
measurable and is in L 2 (dD r ) . 

Next we give an equivalent variational formulation for the system (2.9) and state the 
existence and uniqueness of a weak solution. 

Theorem 2.4. The pair {u,p}, with u G V r and p £ L 2 (Df) is a weak solution for the system 
(2.9) if and only if for every cj) € HQ(D) d such that e(</>) = in D r one has 

N 

/ 2fj,e(u) : e(4>)dx — j pdivd)dx= / / • 4>dx + V^ / fi(x,u) ■ <pds. (2.14) 

JD JD f JD f^ JdD l r 

The proofs of Theorem 2.4 and the following Theorems 2.5 and 2.6 are given in the 
Appendix. 

Theorem 2.5. If {u,p}, with u £ V r and p € L 2 (Df) is a weak solution for the system (2.9) 
then u is the unique minimizer of the functional E r : H^{D) d — > R U {+oo} ; defined by 

N 

E r (v) = / /ie(i>) : e(v)dx - / / • vdx + V] / gi(x,v)ds + I Vr , 

JD JD /=1 JdDl 

where Is represents the indicator function of the set S, defined by 

ifs£S 



Is ^~ \ +oo ifs^S. 

Theorem 2.6. There exists a weak solution {u,p}, with u € V r and p G L 2 (Df) of the system 
(2.9), with u being unique, and p unique up to a constant. 

3 Formulation of the Homogenization Problem 

In this section we formulate the homogenization problem. We focus on a suspension of 
rigid particles in a viscous fluid, periodically distributed, where the size of the particles is of 
the same order as the size of the period. 

Let D be a bounded domain in R with Lipschitz boundary and let Y = I — , — ) be the 
unit cube in Mr. For every e > 0, let N e be the set of all points k G Z such that e(k + 1") 



is strictly included in D and denote by |iV e | the total number of them. Let T be the closure 
of an open connected set with Lipschitz boundary, compactly included in Y . For every e > 
and k £ N e we consider the set T^ CC e(k + Y), where T| = e(k + T). The set T^ represents 
one rigid particle, suspended in the fluid occupying the domain D. 
We now define the following subsets of D: 

Dl= (J T^, D) = D\Dl. 

keN 1 

As in Section 2, D\ represents the union of the rigid suspensions, and D% the part of the 
domain D filled with a viscous incompressible fluid of viscosity [i. Let n be the unit normal 
on the boundary of D% that points out of the fluid region. 

Now we describe the random superficial forces acting on the boundaries of all Tt, for k E N e . 
Let (ft, J-, P) be a probability space and assume there exists a function g : dT x l* 1 x (1 -> 1, 
extended Y-periodically in the first variable to a function defined on Ufcez d (^ + ®T) x R x ft, 
with the following properties: 

(HI) for every fixed z £ M d and for almost every oj E Q, the function s \— > g(s,z,co) is 
measurable with respect to the surface measure on dT, 

(H2) for almost every s € dT and for almost every oj € ft the function z i— > g(s, z, oj) 
is convex and Frechet differentiable at every point z E R , and the Frechet differential is 
continuous on M. d , 

(H3) for every fixed z € M. d and for almost every s € dT the function oj i— > g(s, z, oj) is 
measurable with respect to P, 

(H4) s ^ g{s,0,oj) belongs to L l {dT x 0), 

2 

(H5) there exists 7 G (0, 1] and a function a £ L 2 -t (3T x $1), such that 

Sf(s,«i,w) - g(s,z 2 ,oj) < Ca(s,u)\zi - z 2 \"', (3.1) 

for almost every s € dT, a.s. w € O,, and all 2i,Z2 € Mr. 

On the probability space (£1, T , P), we consider a ci-dimensional dynamical system r, which 
is a family of mappings (Tfc) fcgZ d on $7 that satisfy the following properties: 

(Tl) Group property, i.e. tq is the identity and 

Tfc o n = Tk+i for all k,l £ Z 

(T2) Invariance, i.e. the mappings (?fc) fceZ d are measurable and measure preserving, which 
means that 

P{t^B) = P(B) for all k, I € 1 d and 5 £ 7" 

If the system r satisfies the additional ergodicity property: 

(T3) 

if T k B = B for all k £ Z d then B £ {0, ft} 

then the system is called ergodic. 



For every k G Z , and to any measurable function / in ($7, J 7 , P) we associate the function 
T kf = f ° T k which has the same distribution as /. Therefore, the dynamical system induces 
a (i-parameter group of isometries on L P (Q) for 1 < p < oo. 

A set B G J- is called r-invariant if PfaBAB) = for every k G Z , where A denotes 
symmetric difference. A measurable function / is called r-invariant if for every k G Z d : r^/ = / 
a.s. It can be easily shown that the r-invariant sets form a sub u-algebra of J 7 , denoted by 
I, and a function is r-invariant if and only if it is measurable with respect to I. Then, the 
ergodicity property is equivalent to any of the following: 

B r-invariant implies P(B) G {0, 1}, 

/ r-invariant implies / = constant a.s. 

Given r a ti-dimensional dynamical system, to any function / from L 1 ($7) we associate an 
additive process on finite subsets F from Z with values in L 1 ^): 

S(F,f)(u) = Y,Tkf(u) 

fcGF 

and for every F we denote by A(F, f) the following average of / 

A{F,f)( U ) = -±-S(F,f){ U ). 

\r | 

We now use the pointwise Ergodic Theorem ([25], Ch. 6 or [14] Ch. 8), which states that the 
sequence of averages A(nQ, /) will converge pointwise to the expected value of the function / 
with respect to the a— algebra I, E(/|X) a.s. u G 0,, where Q is a nontrivial cube in Z rf , i.e. 
there exists u, v G Z d with U{ < Vi for every i G {l,...,d} and 

Q = [u,v] = {k £ Z d \ u < k < v}. 

If r is ergodic, then A(nQ, f) converges pointwise to / fdP a.s. w£fl. 

Jn 
We present now an equivalent formulation of the pointwise Ergodic Theorem (see [28] for 
the proof of equivalence) that will be applied throughout this paper. If U G M rf , open, bounded 
with Lipschitz boundary and e > 0, we use the notation 

- = {x G R d | ex G U}. 

Then, if / G L 1 (ri) and r is an ergodic dynamical system acting on Q, the sequence of averages 

A I Z rf n — , / ] will converge pointwise to / fdP a.s. ui G VI. 
V e / Jn 

Define the functions g\ :dT^xR d xn^R by 

g%(s,z,u) = eg(-,z,T k UJ) , (3.2) 

for every k G N e , where e is the size of the periodic structure, the function g satisfies (HI) - 
(H5) and the family (rk)kez d satisfies (Tl) - (T3). We notice now that a.s. u G £1 the function 
(s,z) i->- g^.(s,z,oj) satisfies (Gl) - (G4). 



If f e € L 2 (D) is the body force and if the random superficial forces acting on the boundaries 
of T k are given for every cj £ Cl by 

(/Do, = -V z 5ffc(s,z,r fc u;), (3.3) 

then for every u> £ Cl, the system that we intend to study is the equivalent to the system (2.9) 



- div <r e , 
divu^ 

e«) 

KJ 

cr e ,nds 



in D^, 



f)T! 



/ cr^n x (a; - x e k )ds = \ 

J&Tf J 7 



in D e r , 
on dD e r , 

/ / e cta; + / f k (x,u e ,T k uj)ds, for fc € iV e , 
f € x(x — x e k )dx 
+ / f k (x,u,r k uj) x (x - x e k )ds, for k £ N e , 

JdT? 



(3.4) 



u: 



,dT k 

on 91?. 



where, for every oj € Cl, u^ is the velocity of the fluid, p e u is the pressure and cr^ = 2//e(w^) — 
p^I is the stress tensor associated to u^ and p^. In (3.4), a;| is the center of mass of T k , 
xi = ek + ecco, where Xo is the center of mass of T. 

We define the space V e , equivalent to the space V r defined in (2.11) 

V e = {v € Hl(D) d \ diw = and e(v) = in D e r }, (3.5) 

and we look for a weak solution of the system (3.4), u^ £ V e and p^ € L 2 (D e r). 

We know from Theorem 2.4 that a weak solution of the system (3.4) satisfies the following 
equivalent variational formulation 



D 



2[ie(ul 1 ) : e{cj))dx 



p e u div <pdx 



D% 



r -4>dx+ ^ 



D 



k£N c JdT k 



V u eg[-,u,T k u) ■ (pds 



(3.6) 

for every 4> £ ^o(-^) d such that e(0) = in .DJL We also have, from Theorem 2.6 that the 
solution exists and is unique in V e x L 2 (D%) and from Theorem 2.5 that u% is the minimizer 
of the functional E%, : H^(D) d -> R U {+oo}, defined by 



££( v ) = / (J>e(v) : e(v)dx - / f e -vdx + Y] / e# (-,v,7fcw) ds + /y e (t>). 



(3.7) 



keNe * dT* 



In the next chapter we study the limiting behaviour of the system (3.4) as e — > for 
the suspension of rigid particles in an incompressible viscous flow and determine the effective 
properties of the suspension. We show that the solutions of the system (3.4) converges when 
e —)• 0, to a solution of a system that we will describe later. We take advantage of the 
minimization property of the solution u^, and justify the obtained asymptotics using the T- 
convergence method, introduced by De Giorgi, that implies the convergence of minimizers. 
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4 Homogenization Results 

4.1 Cell Problems 

In order to determine the T- limit of functionals E e , we need first to define the so-called 

CO 

"cell problems" and their solutions. Let S be the subset of R consisting of symmetric 
matrices. For any A S S, consider the following minimization problem: 

min / n[At T - e(v)] : [A1 T - e(v)]dy, (4.1) 

where 

K A = {ve H^ er (Y) d /M | e(v) = A in T, div v = in Y \ T}, (4.2) 

and It represents the characteristic function of the set T, defined by 

. . N f 1 if y € T 

The choice of the space Hp er (Y) /R, which we identify from now on with the subspace of 
Hp er (Y) consisting of zero mean vector fields, is justified by the fact that the functional we 

are trying to minimize v — > / /j,[A1t — e(u)] : [AIt — e.{v)]dy does not change if we add 

JY 

to any v G Hp er (Y) a constant, thus on Hp er (Y) the functional is not coercive and also we 

cannot have uniqueness of the solution. 

A solution for the minimization problem (4.1) thus exists because the functional v —> 

/ /i[Alr — e(t>)] : [A\t — e { v )]dy is convex, continuous in the strong topology of i7 1 (y) rf /M 
JY 

and coercive, and the set Ka is closed and convex. The minimizer xa £ ^A may be also 
characterized by 

J 2ix\A\ T - e(xA)\ : Hxa) - e(v)]dy = 0, 
for every v € Ka, which is equivalent to 

J 2n[Al T - e( X A)} : e(cf>)dy = 0, (4.3) 

for every <fi € Kq. We easily see from here that xa € Ka is unique. 

The formulation (4.3) also implies that xa is the solution of the following Stokes system 
(see also [26] for the formulation of the cell problem) : 

-Axa + Vtia = in Y\T, 

div xa = in Y \ T, 

e( XA ) =AinT, { ^> 

XA eH^ er (Y) d /R, 

with xa & Ka and t\a € L 2 (Y \ T)/R are unique. 

If we extend It and xa by periodicity to the whole Mr, and keep for simplicity the same 
notations, we obtain that xa £ H} (M d ) d as a consequence of periodicity. Moreover, the 
following result holds: 
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Theorem 4.1. For every divergence free vector field G H l {Mr) with bounded support, 

" )|» 



having the property that e{4>)\k+T = 0, V7c G Z d one has 



[ 2u[At T - e{ X A)] : e(cf>)dx = 0. (4.5) 

Proof. Here we follow the idea of [12], Th. 4.28, and cover the support of <p with finitely many 
translated cells of Y, (i^)™ 1; chosen in such a way that d{k + T) and dY{ are disjoint for every 
k G Z d and every 1 < i < m. From the choice of the sets (Yi) 1 £ =1 , it is obvious that for each 
1 < i < m there exists a unique ki G 7L d such that ki + T CC 1^. So, we can associate a smooth 
partition of unity (9i)^i to the covering (Yi) 1 £ =1 with the properties 

di G Cff'P'i), < 9i < 1, 0j = 1 on a neighborhood of fcj + T, Vz : 1 < i < m, 

m 

> ^ 6i = 1 on a neighborhood of supp 0. 



i=l 



If {9i4>)rt represents the vector field defined in Y after we extend Y— periodically 9i<f>, then 
{9icf))^ z is Y— periodic and e((9i<f))^) = on T, so 

„ m ,. 

/ 2n[A\ T - e( X A)} ■ e{(t>)dx = V / 2u[Al T - e(x A )] : e(^0)daj 

m „ 

= Y, M A1 T ~ e{ X A)] ■ e(9i<t>)dx 

i=i ^ 

m „ 

= ^ / 2/z[Al r - e( XA )} : e((9 i( t>)*)dy, 

j=l ^ Y 
and using the formulation (4.4) and the properties of the functions (#j)™ i we obtain that 

„ m „ 

/ 2n[At T - e(x A )} : e(0)da: = V / -r? A div(^) # d<, 
JM d i=1 iy\T 

-n A div(^0) # <% = O. 
i=\ 

D 

We have the following consequence of Theorem 4.1: 

Corollary 4.2. For any A G <S, one /ias %a G H 2 (Y) d and tja G i^ 1 (y \T) and t/ie following 
estimates hold: 

\\xa\\ H 2(y)<i < C\\A\\, \\VA\\m(Y\T) < C\\A\\, 
where the constant C is independent of A G S. 

Proof The first part follows from the regularity theorem ([33], Proposition 2.2), all we have 
to show is that the boundary values of xa on dT and dY are smooth. But the condition 
e {XA) = on T implies that xa G C°°(T) and Theorem 4.1 tells us that outside the sets 
k + T, xa is also infinitely differentiable being the solution of the Stokes system with body 
force. 

To obtain the estimates we notice first the linearity of the map A i— > xa and then the fact 
that A i-> ||XA||.ff 2 (Y) d is a norm on R dxd . □ 
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Now define the bilinear form C on S as follows: 



C[A, B] = J »e( X A) : e( X B)dy. (4.6) 

4.2 Auxiliary Lemmas 

In this subsection we state and prove several auxiliary results needed to show the conver- 
gence results from Subsection 4.3. 

In the r— convergence proof we need first to approximate any vector field from V with 
a sequence of vector fields from V e such that we have also convergence for the functionals. 
Lemmas 4.3 — 4.5 are preliminary results needed to perform the construction. Our strategy 
is to use the solution of the local problems and construct first a sequence in Hq(D) with all 
the required properties to be satisfied up to some small error, and then using a diagonalization 
argument in Lemma 4.5 to obtain the recovery sequence. In Lemma 4.5 we need also to 
estimate how well one can approximate a vector field from H l (D) with one from V € , and 
the result is given in Lemma 4.4. The periodic arrangement of the rigid particles is very 
important here and we take advantage of it using Lemma 4.3 which is an extension type of 
result in the unit cell. 

Lemma 4.3. Let v £ H 1 ^) . Then, there exists a vector field v E Hq(Y) , such that 
e(v) = e(v) in T and \\v\\ H i/ Y \d < C\\e(v)\\ L 2t T yyd. If divv = in T, then v may be also 
chosen so that divi> = inY . 

Proof. Let V\ be the closed subspace of H 1 ^) that consists of vector fields v\ such that 
e(v\) = almost everywhere in T, and V 2 the complement of V\ in H 1 {T) d . On V 2 , we 
consider the norms || • |[i and || • H2 given by: 

1 1 f 2 1 1 1 = ||^2||.ffi(T) d 

and 



IV2II2 = ( / e(v 2 ) ■ e(v 2 )dy\ 



1/2 



The norm | • 1 1 1 is the norm induced and | • 1 1 2 is a norm as a consequence of the definition 
of V2. Obviously we have that | |f 2I I2 < II^Hi for every v 2 G V2, so as a consequence of the 
Open Mapping Theorem, there exists a constant C such that for every v 2 £ V 2 we have 

||V2||l < C 1 1 ""2 1 1 2 - 

Consider v € H l (T) 1 , and let v = v\ + v 2 with Vi G Vi for i = 1,2. Let v 2 be an extension of 
v 2 , v 2 € H^{Y) d , and ||^2||jyifK) d — C||' y 2||i- We define then v to be v 2 , so in T we have 

e(v) = e(v 2 ) = e(v 2 ) = e(«i + v 2 ) = e(v) 

and 

\\v\\ H i( Y )<i < CI 1^2 ||i < C||V2||2 = C ( / e(v) : e(v)dyj 

Now, if div v = in T, let v 2 be the extension obtained as in the previous case, and consider 
the weak solution of the following Stokes system: 

-Aw + Vp = in Y \ T 

div w = div v 2 in Y \ T 
w = on dY U dT 
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1/2 



that satisfies the estimate 

\\ w \\H^(Y\T)d < C|l div ^2||L2 ( y\ T )d. 

We take now v to be V2 — w and v obviously satisfies the three properties. □ 

Lemma 4.4. Let u € #q (D) rf . Then, there exists a vector field u' € V e with the property that 
\\u — u'\\ H ir D \d < C l\\div u\\ L 2t D \ + \\e(u)\\ L 2/ D e\dxd) , where the constant C is independent 
ofe. 

Proof. We first construct u± £ HQ(D) d as the solution for the Stokes system 

-Ami + Vpi =0 in D 

divui = divti in D 
Wi = on dD 

The solution u\ satisfies the estimate 

\\ u l\\Hl{DY < C|l divw llL2(D), 

and now the vector field u — u\ is divergence free and satisfies 

||e(u) - e(«i)|| L 2 (D?) dxd < C [\\e(u)\\ L 2 {D . r) dxd + ||Mi||tfi( D )d 

< C(^||e(w)|| L2(D?)d xd + ||divw|| L 2 (D) J . 

For the divergence free vector field (u — uij^+fiT we a PPly a rescaled version of Lemma 4.3 
and obtain the divergence free vector field 112, fc € H^ek + eY), such that 

e ( u 2,k) = e ( u — u i) in tk + cT 
and 

\\ U 2,k\\L 2 (ek+eY) d + e l I ^ u 2,k I \L 2 (ek+eY) d * d ^ Ce\\e(u - Ui) | | L 2 (efc+eT) dxd , 

where the constant C depends only on T. 

We form the sum of all ti2,fc f° r ah k € N e and obtain tt 2 = > ^2,/b with «2 £ Hq^DY and 

fc€7V £ 
||«2||fl-i(D)d < C||C(« - Wl)|| L 2 (D , )d xd. 

Take u' = u — u\ — 112 and see that 

\\ U - U '\\ul{DY ^ W U ^\\hI{DY + H U 2||j/i(D)d < C(\\dw U\\ L 2 {D) + \\e(u - U^W^^eyxd) 

< C(ll e HllL2(D«)rfxrf + l|divw|| L 2 (D) ), 

and by construction, u' is divergence free and e(u') = almost everywhere on D^, hence 
u' £ V e . U 

Lemma 4.5. Let u £ V and w G O. Assume that for any 5 > i/iere exists a sequence 
u e s £ HQ(D) d weakly convergent to u as e — > swc/i £/ia£ 

limsuplimsup / (divu^) dx = lim sup lim sup / e(lt|) : e(u e s )dx 

= limsuplimsup \E^(u £ s ) — L\ = 0, 

/or some l£t. XTien i/iere exists a sequence u e £ V e , weakly convergent to u as e — > such 
that E^(u e ) — y I as e ->■ 0. 
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Proof. For each vector field u s let u' s be the vector field in V e given by Lemma 4.4 that 
satisfies 

\\ u 8 ~ u'IWh^d)" < C [\\div u e s \\ L 2 {D) + \\e(u e 5 )\\ L 2 {D e )d xdj . 

Condition (4.7) implies that 

lim sup lim sup 1 1 u e g — u'\ \ \ H i r D \ d = 0, 

and using a diagonalization argument ([2], Corollary 1.18) we find a sequence 5(e) — >■ such 
that 

Then, the sequence u e = u'\r t \ is the one we are looking for. □ 

In the second part of the T- convergence result we need to show that the sequence con- 
structed is optimal, i.e. the limit given for any other sequence is greater, and we will use the 
next lemma to prove this. 

Lemma 4.6. Let A € S with zero trace and let w e = Ax — eXA ( _ ) extended to the whole D 
by periodicity. Then, for every sequence u e G V e : u e — *■ ase-> 0, one has 

lim sup / 4>e(w (L ) : e(u e )dx = 0, 
e^o Jd 

where <\> £ W 1,co (D) compactly supported in D. 
Proof. We have that 

e(w e ) : e(^) = M™! ■ e(u e ) + £ e tJ (w<) (^-u) + ^-u{\ , 



then 



lim sup / 4>e(w e ) : e(u e )dx = lim sup / e(w e ) : e(cpu e )dx. 
e->o Ju e-5>0 Ju 



Hj 
The sequence 4>u e is weakly convergent to and 



||div0M e || L2(D) < C||V0|| L ^ {D)d ||w e || L 2 (D) d, 

/ e((j)u e ) : e((f>u e )dx < C\\V<f>\\ 2 Loo{D)d \\u e \\ 2 L2{D)d . 

For every <j)u e we apply Lemma 4.4 and construct v e G V e such that 

\\4>u e - v e \\ H i {D)d < c(||divw|| L 2 (D) + \\e(u)\\ L 2 {D e)dxd) < C\\V(f)\\ LX{D)d \\u e \\ L 2 {D)d , 

with a constant C independent of e. We notice that due to Theorem 4.1 

/ e(w e ) : e(v € )dx = 0, 
Jd 

and using the strong convergences to of 4>u e — v e we obtain that 

lim sup / (f>e(w e ) : e(u e )dx = lim sup / e(w e ) : e((pu e )dx = 0. 
e^o Jd e-^o Jd 



□ 
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We derive in the next lemma an important estimate satisfied by the functions gl used in 
the definition (3.3) of the random superficial forces. 

Lemma 4.7. Let the functions g e k : dT k x R x Q, — > R be given by (3.2). Then for every 
u, v £ H l (D) d the following estimate holds: 

/ g e k (x,u,uj)ds- / gl(x,v,uj)ds<C e {u)(\\u-v\\l 2( )d + e 1 \\Du-Dv\\l 2( )dxd ) , 
JdD% JdD% v K ' K ' ' 

(4.8) 
where there exists a constant C such that C e (u) — > C a.s. uj £ Q as e — > and 7 is the Holder 
constant given by (3.1). 

Proof. We use first the property (3.1) and notice that one can assume the function a is non- 
negative. We obtain that a.s. w E £1 we have 



/ gl{ x ^ u ^) ds -j 9k( x ' v i UJ )ds= Y] / e(g(-,u,T k uj) - g(-,v,T k uj))ds 
JdDi, JdD% k TN- JdT t e e 

> ^ / ea i — , rjtw) |it — t>| 7 ds. 



Then from Holder's inequality and the change of variables — = s we obtain 

e 

2 — "/ 7 

a ( — , TfcCj) lit — vPds < I / a ( — ,rjtw ) 7 (is / Itt — t>| 2 ds 



9T| v e 7 W9T| v e 



7 
^ / \ _ 

(n-l)(2- 7 ) / /' 2 \ 2 



2=7 / . \ 2 



e' 2 ( / a(s,Tfew) 2 -T (is ) j / |« — t>| 2 <is 

which implies that 

/ g%{x,u,u)ds- / fl£(x,v,a;)ds < V / ea ( -,r fc w) |« - "u| 7 ds 
J9DJ; JdD% ,~Z T . JdTf v e 

-. (n-l)(2- 7 ) ^-^ / /" _2_ V 

<e 1+ 2 \ / a (*,7iw)»-Tds 

„ „ 1 I ("~l)(2-7) / v-^ f / s-?- , 

< Ce 1+ 2 V / a(s,r fe w)2-7 d s 

We obtain in the end that 

eg 1-,u,oj) ds — / eg ( -,v,u>) ds < 

2-7 

\*t^ |7V 1 JdT J \k^ Jd n , 

2 / f \ 2 



fcG7V £ ' 


lffT£ 


V t 








2-7 / 




— v\ 


2 ds 


2 

1* 




) 




f 


\u - 


v\ 


2 ds 



<Ce2 V— — / a{s,T k uj)^ ds) V / \u - v\ 2 ds 

<CeiC e (w)( V / \u-v\ 2 ds) 



vfceA^ U1 k 
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where, by the Ergodic Theorem C e (io) — > \\a\\ 2 as e — > a.s. we!!, and from the 

property (H5), Hall 2 is finite. 

Applying a rescaled Poincare inequality 



|tt — i)| 2 iis < e 1 [ I \u — i>| 2 da: + e 2 / |Z>u — Di>| 2 d:K 
one has 



/ eg (-,u,Ld) ds— eg ( -,v,u>) ds < C e (oj) [ / \u — v\ 2 dx + e 2 \Du — Dv\ 2 dx\ , 

which yield the desired estimate. 

In order to prove our main convergence result we need an improved version of the pointwise 
Ergodic Theorem that we state and prove in the following theorem. This theorem can also be 
viewed as an extension of the pointwise Ergodic Theorem proved in [28] with weights of the 

form j- u{x)dx added in computing the averages, where u G L 1 (C7), and U C M. d is 

Jek+e[0,l] n 

open, bounded with Lipschitz boundary. □ 

Theorem 4.8. Assume 1 < p < +00 and h : M. m x Q, — y R a function such that oj 1— > h(0,co) 
belongs to L l (Q) and for every Z\,z 2 G M m and a.s. ui G £1 



h(zi,co) — h{z2,oj) < \z\ — z s \ p a(io), 



for a function a G L p (fl). 



Then, on a set of full probability the sequence e n > h I + u(x)dx,TkUJ con 






feeZ «n^ v^[o,il" 



verges to / h(u(x),uo)dPdx for any U C K n open, bounded with Lipschitz boundary and any 



ue L l {U) m . 

Proof. Obviously, for every z G K m , the function u 1— > h(z,uj) belongs to L 1 ($7). Let Q,' 
be a set of full probability such that for every ui G Q' and for any z G Q m , the sequence 

y e n h(z, TkUj) converges to / / h(z,uj)dPdx, for every U G R n , bounded with Lipschitz 

, „ ,, Ju Jn 

k£Z n nf 

boundary. Moreover, 

^2 e n h(z 1 ,T k u!) - ^2 £ n h(z 2 ,T k u>) < \z\ - z 2 \~p~ ^2 e " a ( r fe^)i 

therefore passing to the limit e n \. h ( z > T k w ) dx converges to / / h(z,u))dPdx for any 

, „ v Ju Jn 

U C R n open, bounded with Lipschitz boundary, any z G R m , and a.s. u> E CI. 

Assuming now that u G C(U) m , where U is as in the assumptions of the theorem, and let 
xq G M. n be such that 

/ / h(u,uj)dPdx = / / h(u(xo),uj)dPdx. 
Ju Jn Ju Jn 
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Then, 



lim sup 



2_] e n h(u(ek),TkOj) — / / h(u,cj)dPdx 
"~~ ~" Ju Jn 



fcez™n^ 
a.s.. The continuity of u implies that 



p-i 



< l|o||x,im)(supw-infw) * C n (U), 



u u 



lim V^ e n h(u(ek),TkOj) = h(u,oj)dPdx, 

e 

a.s. u> £ Q. In order to prove the theorem, we use the separability of L l (U) m and the density 
of C(U) m in L l (U) m . We first notice that we can show in the same way that 



lim 



2_] £ n h j- udx,Tk^> 1 = / / h(u,u>)dPdx, 

:»„£ \Jek+e[0,l] n J JuJa 



k&z n nf 
for u € C(U) m . If u € L^U) 171 and v € C(U) 



2_] £ n h t udx, TkCJ — 2_] e n h + vdx, T k 0J 



< 



p-i 



V] e n If (u - v)dx a{T k ixj) < 



feez™ni 



np 



kez n n 
This implies that a.s. to & Q 



2_] / (u — v)dx 



^2 a ( T k UJ ) 1 



\ 



,fcez»n- 



lim sup 






feeZ"ni 



e n h 



{£ 



fc+e[0,l] n 



udx, T^OJ 




u Jn 



h(v,u>)dPdx 



<C||u-v||,? 



L i (uy 



Now we let v converge to u in L (U) m and obtain the result. 



□ 



We use the Ergodic result given by Theorem 4.8 and the estimate given by Lemma 4.7 to 
derive the following two limits, that essentially show how the random, highly oscillating forces 
acting on the boundaries of the rigid particles, become through homogenization, deterministic 
forces in the volume. 

Theorem 4.9. Let u and v be vector fields in H l (D) d and let u e and v e be two sequences 
such that u e —*■ u and u e —*■ v in ff 1 (L>) as e — > 0. Then, on a set of full probability we have: 

i ) lim y / eg I - , u e , r k uj ) ds = / / / g (s,u,co) dPdsdx. 
e-s-O.^ .Ifrre Ve / JoJ&rJn 



fcgATe J 9Tk 






ii) lim , 

k£N f - al k 



tV u g (-,u e ,TkOj) ■ v € ds = / / / V u g(s,u,oj) ■ vdPdsdx. 

e ' J D JdT Jn 
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Proof, i) According to Lemma 4.7 



lim sup 



/ eg (-,u e ,T k uj) ds- V] / eg ( -, u,T k u; ) ds 



k£Nl Jdl k 



fceAfE ^ui k 



e^O 



< limsupC e (w) ( \\u e — u\V L 2t D }d + e^lVit 6 — V« 



lL 2 (D) a 



(4.9) 



e-K) 

We will show now that 



< C lim sup | \u e s — u\ \l 2 , D)d = 0, a.s. ueH. 



B = lim sup 



y / tg ( -j it, r k ui\ (is — / / / g{s,u,u))dPdsdx 
kTN^ dT k e JDJdTJn 

k€N eJO-l k 

y e / ^f (/c + s, w(e/c + es),Tfcw) (is = N^ e / 5 (s,u(ek + es),T k u) ds. 



0. (4.10) 



hew 
Denote by u\ = j- 

Jek+tY 



fcGiV' 



udx. then 



g(s,u(ek + es),T k uj)ds - g (-,u k ,T k uj) ds < C(e) I a(s,T k uj)ds, 
dT JdT JdT 

where C(e) — > as e — > 0. Using the Ergodic Theorem, we obtain that a.s. w G J7 



lim 



^ e d / g(s,u(ek + es),T k uj)ds- ^ e d g (-,u e kl T k uj) ds 

keN€ JdT keNe JdT 



So 



i? = lim sup 



: lim sup 



V / e d g(-,u%,T k uj)ds- / g(s,u,u)dPdxds 

k£Nc JdT JD JdT Jn 

z2 e<1 9 ( T udx,r k uj ) — / g(u,oj)dPdx 



where g is defined on R x Q, by g(z,uj) = / g(s, z,uj)ds. From Theorem 4.8 we obtain 

JdT 

that B = 0. 

ii) Using the definition of the subdifferential we have that for every t £ R 

£V U # (-,w e ,T fc u;J -v € < g (-,u e + tv e ,T k u}\ - g (-,u e ,r k u} 

therefore 

t^2 / eV ^ (-,w e ,TfcwJ • v e ds < ^2 / (eg y-,u e + tv e ,T k u?j -eg \-,u e ,T k u) J c/.s. 

keNe JdTl '' < ■ keNf _ JdT^ e 
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Using the first part, we obtain that 

limsupt y / e'Vug i-,u e ,TkUi) ■ v e ds < / / / (g (s,u + tv,uj) — g (s,u,u>)) dPdsdx 

< / / / ~V u g (s,u + tv,uj) ■ tvdPdsdx. 
Jd JdT Jn 

We take first t > 0, and then £ < and let it converge to 0, and then apply the Dominated 
Convergence Theorem and obtain that 

lim > / e'Vug (-,u e ,TkUj) ■ v e ds = lira / / / ~V u g (s,u + tv,u>) ■ vdPdsdx 
e ^°k7N^ dT k ^JDJ&rJck 

= I 'Vud (s,u,u>) ■ vdPdsdx. 

Jd JdT Jn 

□ 

4.3 Convergence Results 

Assume there exists / E L 2 (D) d such that / e — *■ / and introduce now the functional 
E* : H 1 ^)- ->RU {+00} defined by 

E*(v) = / (pe(v) : e(v)+C[e(v),e(v)])dx- fvdx+ / / g(s,v,u)dPdsdx+I v (v). 
Jd Jd Jd Jbt Jn 

(4-11) 
Properties (H4) and (H5) of the function g imply that there exists a constant C such that 

+ C\z[< 

(412) 



g(s, z,u>)dPds 
dT Jn 



g(s, 0,cu)dPds 
or Jn 



< 

<C + C\zp, 



for every z € Mr. Thus E* is well defined on H l (D) as a consequence of the property (H2), 
is strictly convex and takes finite values on V. Moreover (4.12) implies that there exist two 
constants A > 0, and C such that for all v £ H l (D) we have 

E *(v) >-C- C\\v\\l 2{D)d - C\\v\\ LHD)d + A\\v\\ 2 L2{D)d , 

therefore there exists an unique minimizer u* € V for the functional E* . It can be shown 
similarly as in Theorem 2.5 that the minimizer u* will also satisfy the following equivalent 
weak formulation: 

/ 2fie(u*) : e(<j>)dx + [ 2C[e(u*),e(<f>)]dx = [ f ■ 4>dx + f /*(«*) • <f>dx, (4.13) 
Jd Jd Jd Jd 

for every </> € V, where we denoted by /* the vector field on Mr defined by 

f*(z) = -[ [ V z g(s,z,Lu)dPds. (4.14) 

JdT Jn 

The linear functional defined on H l (D) d 

<t>\-+ I 2fie(u*) : e(<p)dx + / 2C[e(u*), e(<p)]dx - [ f ■ <j>dx - [ f*{u*) ■ <$>dx 
Jd Jd Jd Jd 
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being on V, implies the existence of p* € L 2 (D) such that 

f 2fie(u*) :e(<f>)dx+ [ 2C[e(u*),e(4>)]dx - I p* div cj>dx = [ f ■ d>dx + I f*(u*)-d>dx, 
Jd Jd Jd Jd Jd 

(4.15) 

for every <f> G H^{D) d . 

The pair {u*,p*}, u* € H 1 (D) d , p* € L 2 (D), satisfying equation (4.15), is a solution for a 

Stokes system that we describe now. Given the bilinear form C, we denote by /x* A the unique 

element from S such that 

fiA: B + C[A,B} = n*A: B, (4.16) 

for all B £ S. We also denote for any u € H l {D) and any p € L 2 (D) by cr*(u,p) the stress 
tensor 

a*(u,p) = 2^*e(u) - pi. (4.17) 

Equation (4.13) can be written equivalently in the following way: {u*,p*} is the solution for 
the Stokes system 

-div (T*(u,p) = / + /*(«) in A 

divw =0 inD, (4.18) 

u = on 3D. 

We will show next that E* is the limit functional for the sequence of functionals (E£j) e >o 
and u* , the unique minimizer for E* is the limit solution, i.e. E* is the T- limit for the 
sequence (E^) e> Q and u e w —*■ u* in V, a.s. ui € fi. Then, we have the following result: 

Theorem 4.10. The sequence of functionals (E^) e> Q T- converges to the functional E* in the 
weak topology of H 1 (D) d a.s. uj € ft, where E* : i7 1 (D) d -> R U {+00} is defined by (4. 11). 

Proof, i) For the first part we need to show that 

r-limsup££ <E*, 



a.s. uj e ft. 



We take first u € Ffl W 2 '°°(D) d and let 5 > 0. Then, there exists a partition of D, (A)il\ 
such that for all 1 < i < N$, Di is Lipschitz and 

essosce(ti) < 5, 
A 

where essosce(w) denotes the essential oscillation of e(tt) over the set Dj, i.e., the smallest 

Di 

number t > such that ||e(u(aji)) — e( , u(a?2))[| < t, for almost all xi and X2 in Z?j. 

We define Aj = + e(u)dx and let wf(a;) = eXAi ( — ) extended by periodicity to the 

whole D{. For every 1 < i < Ng we denote by Di^ the set of points from Di with the distance 
to the boundary of Di at least 5, and let fag be a smooth cutoff function supported in Dj, 

C 

such that < &,$ < 1 in Di, equal to 1 in Di t s and such that \\^(j)i,8\\L°°{Di) d — ~T- Let us 

define w| by 



u\=\u 



y <\>i^u\ — ^ w as e — > 0. 
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We now show that the sequence u\ satisfies properties (4.7) a.s. w € f2. It is easy to check 
that || divit^| \l 2 (d) -^ as e ^ 0. Also 



e (Ux) : e [us 



< 



)dx = Y] / e(u|) : e(u € 5 )dx 

i=1 Jd^d, 

V 2 / e(w - w|) : e(w - itf)cfcc + V 2 / e(w t e (l - ry^)) : e(w, e (l - (f> ijS ))dx 

r Ns r N 5 

<C 5 2 dx + V C / e«) : c(«|)(l - <M 2 ^ + 12° \ u i\ 2 \^M 2dx 

JD ~!y Jd%C\D % ~[ Jd { 

KC^ + f^cf e(uf):e(uf)dx + J2%2 [ \<\ 2 dx, 

i=1 JDi\D ijS i=1 d JD, 



which implies that 



Finally, 



lim sup lim sup / e(w|) : e(u e s )dx = 0. 



K( u s)- E *( u )= / A»c(«|) : e(«|)dx - / f € ■ u%dx + V / e^ (-,u§,Tfcw) cte 

Jd Jd cr? T JdTi v e / 



fcew u± k 



l (fj,e(u) : e(w) + C[e(u) , e(u)]) dx + / / • udx — / / / g(s,u,cj)dPdsdx, 
Jd Jd JoJdTJn 



hencefore 



where 



lis = h m SU P 



limsup |.E£(u|) - £?*(u)| < h, S + h,8 + h,S, 



jue(u|) : e(u € x)dx — / (//e(w) : e(w) +C[e(w),e(w)]) da; 



D 



I26 = hm sup 

e->0 



/ / e • u e gdx — f • udx 
Jd Jd 



and 



Izj = lim sup 



^ / fry I -.u};.n-,c 
keNe JdT* 



ds — / g(s,u,u>)dPdxds 

Jd Jar Jn 



It is easy to see that I2S = 0, and according to the first part of Theorem 4.9 I3 ,5 = as 
well. 

To estimate Ii^ we make the following elementary calculations: 

/ jtxe(u|) : e(u e s )dx — / (pe(u) : e(tt) + C[e(w), e(tt)]) da; = 
jd Jd 

2_)\ / l^e(u — (f>i xu\) : e(u — (f>i xu\)dx — j (/j,e(u) : e(u) + C[e(u) , e(u)]) dx ) = 
~{ \JDi Jd, J 
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y~] ( / -2fj,e(u) : e{^ s ul)dx + / ne(4>^ s u\) : e(<j)i jS u e i )dx - C[e(u),e(u)]dx ) . 

i=1 \JDi JD l JD t J 



Since 4>i^u\ — *■ as e — > 0, 



i^i^uf) : e(4)iju\)dx — I C[e(u), e(u)]dx 



o, 



I\ s < V^ lim sup / //e 

i=1 «->0 Jd, 

NJlimsup / ^i^e{u\) : 4>i^e{u\)dx — / C[e(w),e(w)]rfx 

r-r' e-+0 Jd, ./D; 

lim sup C / e i u i) '• e.{u\)dx + lim sup 

e->0 JDi\D is ' e^O 



< 



< 



(jue(iif) : e(«|) — C[e(w), e(w)]) da; 



D; 



£ 

1=1 
From the construction of wf, / e(uf) : e(uf)dx converges to / C[e(% J 4.), e(% J 4 i )](ia; and 

C[e(x^), eixA)} - C[e(«), e(«)] = C[e(«) - e(x^), e(«) + efoO] < C<5, 

almost everywhere in Di, so 

^ - 

h 6 < C yj lim SU P / e(w| ) : e{u\)dx + C<5. 

^ e->0 JD,\D itS 

We obtained for u G V n W^ 2,oc (Z)) and <5 > the sequence u$ weakly convergent to u in 
Hl(D) d that satisfies properties (4.7), i.e. 

lim sup lim sup / (divtt|$) dx = lim sup lim sup / e(ii|) : e(u e s )dx 

6^0 e^O JD S^O e^O J D% 

= limsuplimsup \E^(ug) — E*(u)\ = 0, a.s. uj € Q.. 

To this sequence we apply Lemma 4.5 for L = E*(u) and obtain a sequence u e € V e , weakly 
convergent to u such that E^(u e ) —> E*(u) a.s. u) € Q. 

We notice that we have proved so far that 

T - lim sup El (u) < E* (u) , 

in the space V equipped with the sequentially weak topology of HQ(D) d , for any u in a 
dense subset of V with respect to the strong topology of H(j(D) d . We now use the fact that 
r — lim sup E^ is lower semicontinuous by the definition, and that E is continuous. 

Take u € V, and let (v,k)k>i be a sequence strongly converging to u in H^{D) d , such that 



r-limsup£E(ti fc ) <£*(«*) 



Setting fc^oowe obtain 



r — lim sup E^ (u) < lim inf T — lim sup E^ [uk ) < lim inf E*(uk), 
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which implies that 

r-limsup^(w) <E*(u), 

for every u G V . 

ii) In the second part we prove 

T - liminf Ef, > E, a.s. wefl. 

For this, we have to show that for every u € H l (D) d and for every sequence u e weakly 
convergent to u we have 

E*(u) < liminf Ef,(u e ). 

Take first w a smooth vector field, w £ V H W 2 '°°(D) and for a small positive 5 > select a 
partition (Di)i = \ ) N 5 of D, such that for all 1 < i < iV$, Dj is Lipschitz and 

essosce(w) < 8. 

As before, Ai = j- e(w)dx and let wf(x) = A{X — exA t I — ) extended by periodicity to 

JDi V e / 

the whole -Dj. Then w| — *■ wi as e — > 0, where Wi(x) = A{X. Take also 0j a smooth cutoff 
function supported in Di with values in between and 1. We have: 

e(u e (x)) : e(u e (x)) > e(wf(x)) : e(w\(x)) - 2e(w e i (x)) : e(u e (x) - w\(x)) 

for almost all x £ Dj. We multiply [Mpi and integrate over Dj and obtain 

/ fie(u € ) : e(u e )dx > / [Mf>ie(wj) : e(wf)dx — 2 / fi<f>ie(w':) : e(u e — w\)dx. 
J D % J D % JDi 

From the choice of Wi we have that 

lim / /MJ)ie(wi) : e(w^)dx = / n4>iAi : Aidx + / <f>iC[Ai, Ai]dx, 

e ^ JDi JD % JDi 



thus 



liminf / fie(u € ) : e(u e )dx > / [ufiiAi : A^dx + / <piC[A{,Ai]dx 

e ^>° JDi J D % JDi 

— 2// lim sup / 4>ie(w' i i ) : e(u € — w\)dx 



Assume for now that 



lim sup / (pie(Wi) : e(u e — wf)dx 

e->0 JDi 



is bounded from above by C||e(it — Wi)\\ L 2t D .\dxd, where the constant C is independent of Di 
and fa. We now sum over all 1 < i < Ng which yields 

liminf/ fie(u e ) : eiu^dx > > / (n^Ai : Aj + faClAi, AA) dx — C\\e(u) — Ai\\ L 2( D .\dxd. 
e ^° Jd T^i JDi y %) 
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We let now fa to converge increasingly to one in each Di apply the Lebesgue Beppo-Levy 
Theorem to obtain 

liminf / /ie(w e ) : e(u e )dx > V^ / (uAi : Ai + C[Ai,Ai]) dx — C\\e(u) — Ai\\ L 2( D .)dxd. 
^° Jd ' ^JDi 

Set S — > and the uniform continuity of e(w) in D implies that 

limV( / fiAi : Aidx + / CLAj, AAdx ) = / fie(w) : e(w)dx + / C[e(w),e(w)]dx, 
s ^° ~\ V A JDi J Jd Jd 

and 

N s 

]im^ ||e(w) - Ai\\ L 2 {Di yxd = \\e(u) - w\\ L 2 {D yy.d. 

Since u e —>■ u as e — > we have according to the first part from Theorem 4.9 that a.s. u£fi 

liminf E^(u e ) > E*(w) — C\\u — w\\ H i (D)d . 

Now we make w converge strongly in H l {D) to u and obtain the desired inequality. 
It remains to show that 

limsup / faeiwf) : e(u e — w\)dx < C\\e(u — Wi)\\ L 2( D .\dxd. 
e^o Jd, 

For every symmetric matrix A with zero trace, and every <fi £ W l, °°{D) compactly supported 
in D, define the functional g:y->RU {+°°} by 



q(v) = sup limsup / <fie(w € ) : e{v e )dx, 
v^ev e^o Jd 



where w e (x) = Ax — ex ( — ) • It is clear from the definition that q is sub linear, i.e. 

q(Xv) = Xq(v) for every positive A and every v € V, 

and 

q(vi + v 2 ) < q(vi) + q(v 2 ) for every vi, v 2 € V. 

It is enough to show that q is continuous with respect to the strong topology on H 1 (D) , i.e. 
there exists a constant C such that for every v G V 

q(v) <C\\v\\ H i {D)d , 

and the constant C is proportional with ||0||l°c(d) an d with ||A||. 

But q(v± + v 2 ) < g(t>i) + limsup / <pe(w e ) : e{v\)dx : for every V\ , v 2 € V, and any weakly 

e^o Jd 
convergent sequence v\ to v 2 . But from Lemma 4.6 q(0) = 0, and in the first part we showed 

that for any v € V, there exists a sequence v e — *■ t> as e — > such that 

limsup ||e(V)|[ L 2 (D) dxd < C||«|| H i (D)( j, 

for a constant independent of v, which means that for every v £ V we have 



g(v) < limsupC||0e(t« e )|| I/ 2(md X d||'y|| fl -i (D) d < C| |0| |z,=«(z3) | |v| | jH -i ( _ D) rf 1 1>1| |. 

e-»0 



D 
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The following corollary provides us with the existence of the homogenized solution, and 
the corresponding energy convergence. 

Corollary 4.11. Almost surely u> G 0, the sequence u^ —>■ u* in V and E^u^) — > E*(u*) as 
e — > 0, where u* is the unique minimizer of the functional E* . 

Proof. We only need to show that a.s. uj £ Vt the functionals ££, are equicoercive. After 
that, from the Y— convergence of functionals E^ to E* and equicoercivity of E^ we deduce 
the convergence of energies as well as the convergence of minimizers ([13] Theorem 7.8 and 
Corollary 7.17), therefore a.s. uj € f2 we obtain that u^ — *■ u* in V and E e (ulj) — > E*(u*) as 

To show equicoercivity we use first (4.10) and obtain a.s. uj € fi the following estimate for 
e small enough and a constant C independent of e: 



E 

k&N* 



eg (-,v,T k u)) ds 



<C+ / / g(s,v,oj)dPdxds, 

Jd JdT Jn 



for any v 6 i? 1 (D) c( . Then we use (4.12) and derive 



E 



k£N e " &TI 



f.9 



,V,T k UJ) ds 



<C + 

'D JdT Jfl 

<C + C\v\~<, 



g(s, 0,uj)dPdxds 



+ r,\\i,\V 

^ u ll t; llL2(D)a 



for a constant C independent of e. Finally, the definition (3.7) of E^ gives us that a.s. uj £ f2, 
there exist two constants A > and C, both independent of e, such that for all v € H 1 (D) d 
and for e small enough we have 



E e u (v) >-C- C\\v\\l 2{D)d - C|M| L 2 (D)d + A\\v\\ 2 L 2 {D)d , 
which implies the equicoercivity of E^. 



□ 



4.4 Corrector Result 

The V- convergence result provided by Theorem 4.10 together with the Corollary 4.11 

gives us only weak convergence in H l (D) of the solution u^ to the homogenized solution u*. 
Next, we formulate the first order corrector result to improve the weak convergence to strong 
convergence. 

First, we define for every 1 < i,j • < d the matrix Eij € S componentwise by 



1 



(Eij) hi = 7T (Sikfyl + &ilb~jk) , 



(4.19) 



for every 1 < k, I < d and redenote for simplicity Xij '■= Xe 13 - Since Xij £ ^-e^ we use the 
definition (4.2) of the sets Ka, and we obtain that e(xEi) = Eij in T. The definition (4.19) 
will imply that 

divXij = SijtT- (4.20) 

Theorem 4.12. Assume that the homogenized velocity u* € V fl W 2,oc (D) d , then a.s. uj € ft, 

( d 

the sequence u^ — \ u* — \, e Xij (~) &ij( u *) — > in H l (D) d as e — >■ 0, where Xij lri 

extended Y - 'periodically to the whole D. 
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Proof. We notice that we need the assumption that u* G V l~l W^ 2 '°°(.D) in order for r e = 

d 

y^ exjj (-) eij(u*) G H 1 (D) d . First, we see that 

and because ti* G VF 2,00 (D) d we have 

d 
IK, - («* - r e )\\ L 2 {D)d < IK - u*\\ L 2 (D)d + e||«*|| w - 2 ,oo( r) )d Yj IIx« (-) lli 2 (D)d 

which implies that m^ — (u* — r e ) — > in L 2 (D) d as e — > 0. If we are able to show that 
e(ulj — (u* — r e )) — > in L 2 (D) as e — > 0, then applying Korn's inequality we obtain the 
desired result. 

To show this we use results of Theorem 4.10. An elementary calculation gives 

e«-w* + r e ) : e«-u*+r e ) = e«) : e«) + e(-u*-r e ) : e(M*-r e )-2e«) : e(u*-r e ). 

The variational formulation (3.6) yields 

/ 2fie(u e J : e(u e Jdx + V / eV u g ( -,u e u ,Tku) ■ u e u ds = / f e -u e u dx, 

Jd keN< JdT * Ve J Jd 

then, using the second part of Theorem 4.9 

lim / 2^e(ti^) : efa^dx = / fu*dx — / / / *V u g (s,u* ,oj) ■ u*dPdsdx 

= I f ■ u*dx + I f* («*) • u*dx 
Jd Jd 

Also, because r e —*■ in H l (D) d as e — > 0, 

lim / 2fie(u* — r e ) : e{u* — r e )dx = / 2fie(u*) : e{u*)dx + lim / 2fie(r € ) : e(r e )dx 
€ ~>°Jd Jd ^°Jd 

But 

lim / 2/ie \exij (-) eij{u*)j : e (eXkl (-) e k i(u*)j dx = 

lim/ 2/ie(xij)(-) ■■ e(xki) (-) eij(u*)e k i(u*)dx = / 2C[E ij ,E k i]e ij (u*)e k i(u*)dx, 
e-¥Qj D Ve/ Ve/ Jd 

hence 

lim / 2fie(r e ) : e(r e )dx = V / 2C[E ij ,E k i]e ij (u*)e k i(u*)dx = / 2C[c(«*), e(u*)]cte 
^°./d ij^i^ " ^ D 

which implies that 

lim f 2 / ue(w* - r e ) : e(u* - r e )dx = [ 2/*e(u*) : e(w*)da; + / 2C[e(«*), e(u*)]da 
£ ->°./d Jd Jd 

= I f ■ u*dx + / /*(«*) • w*dcc, 
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/a; 



after we use the variational formulation (4.13) for u*. Let us show now 

lim / 2fie(ul) : e(u* - r e )dx = f f ■ u*dx + [ f*{u*) ■ u*dx. 
^° Jd Jd Jd 

This is equivalent to showing that 



lim sup lim sup 

8^0 e^O 



I 2^e«) : e(u* - r e )dx - [ f ■ u*dx - [ /*(«*) • u*dx 
Jd ~ Jd Jd 







because the sequence (it^) does not depend on 8. For a fixed 8 > we denote by D$ the 

set of points from D with the distance to the boundary of D at least 8 and select a smooth 

cutoff function (f>$ supported in D, such that < 0,5 < 1 in D, equal to 1 in Dg and such that 

C 
\\^4>8\\L°°(D) d < ~f an d see that 

/ 2jxe{ulj) : e(u* — r e )dx = / 2/ie(u^) : e(u* — r £ )<psdx+ / 2/ze(tt^) : e(u* — r e ){\ — <j)$)dx. 
Jd Jd Jd 

We estimate the second term, and a.s. oj € Q 



I 2 F «):e(ii*-r ( )(l-^ 

JD 



< CWu^Wui^yWu* - r e \\ H i( D \ Ds y 



< ( ' ! \\ u *\\imD\D S y + 2^ \\ e Xij y-J llffi(D\D 5 )rf 



»,j=i 



which implies that 



lim sup 



/ 2 M e«) : e( 



u*-r e )(l-<^)cta: 



< C ! llw*|| H i( D \D 4 ) 






Therefore 



lim sup lim sup 



/ 2^e«) :e(«*-r e )(l-^)daj 



It is sufficient to show that 



lim sup lim sup 

5^0 e->0 



/ 2^ie{ulj) : e(u* — r e )(f>gdx — / fu*dx— / f*(u*)-u*dx 
Jd j ~ Jd Jd 

We perform the following calculation: 

/ 2^e{ulj) : e(u* — r e )(psdx = / 2^ie{ulj) : e(w* — r e <f>s)dx+ 
Jd j ' Jd 

J 2^e«) : e(u*)(<Ps - l)dx + j^Y. e ^O (^ r i + f^<) dx > 



0. 



and estimate the second term 



lim sup 



2^e«) : e(t**)(& - 1) 



D 



< C|| u *llj/ 1 (D\D 4 ) d - 
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Because r e — > in L 2 (D) d ase-)-0 



lim sup 



~ d 

/ 2/J J2 e iJ« 
JD i,7=l 



— — r H — r- 

5xj 3 dxj % 



0. 



It remains to show that 



lim sup lim sup 

5^0 e^O 



/ 2/j,e{ul):e{u*-r € fa)dx- / fu*dx- / f*(u*)-u*dx 

JD '" JD JD 



0. 



For the sequence u* — r e (ps we calculate || div(ii* — t ,£ </'j)||l 2 (d) an d ||e(w* — ^ £ 0<5)||L 2 (D^) fl 
div(w* — r 6 ^) = — (frsdivr 6 — r e ■ V</>5 



u 



= -fae Y^ Xij (-) • Veij( 

According to the property (4.20) in any set of the form ek + eT, for k € Z' 
d d 

X^ divx^ f -J ey(it*) = ^ 5ijeij(u*) = divu* = 



X] div Xy (-) e *i( w *) - »* £ • V( i 

d 



«J=1 



*J=1 



and because div^y = outside the sets ek + eT, we obtain that \, divxij ( - ) &ij{u*) = 
in D. Thus 



lim || div(w* - r e ^)|| L 2 ( m = 0. 



Also in Di 



e(u* - r € fa) ■ e(u* - r e fa) = e(u* - r e ) : e(u* - r e ) + e(r e (l - fa)) : e(r e (l - fa)) 

< e(u* - r e ) : e(u* - r e ) + Ce(r e ) : e(r e )(l - fa) 2 + C|r e | 2 |V^| 2 

< Ce 2 + Ce(r e ) : e(r e )(l - fa) 2 + C|r e | 2 |Vefo| 2 . 



Hence, 



lim sup lim sup / e(u* — r e (f>s) : e(u* — r e (f>$)dx < lim sup / e(r e ) : e(r e )dx = 0. 

5^>0 e^0 JD% 5-+0 JD\D S 

We apply Lemma 4.4 to the sequence (u* — r e (j)s) and find a sequence u^ € V e such that 



This means that 

lim sup lim sup 

5^0 e^0 



lim sup lim sup 



lim sup lim sup 



limsuplimsup ||w* — r e <ps — u\\\ H ii D \d = 0. 

5^0 e^0 



/ 2/j,e(ul) : e(u* - r e fa)dx - / fu*dx- / f*(u*)-u*dx 

JD " JD JD 

/ 2fie(u e J : e(u e 5 )dx - / fu*dx- / f*(u*)-u*dx 

JD ~ JD JD 

[ f e ■ u%dx -If- u*dx + 

JD JD 



< 



lim sup lim sup 

5^0 e->0 



after application of the second part of Theorem 4.9. 



□ 
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5 Conclusions 

The current paper provides a homogenization-type result stated in Theorem 4.10 and 
its corollary for a suspension of rigid particles in the presence of highly oscillatory random 
interfacial forces. In particular, it states that the solution u^ to the family of problems (3.4) 
corresponding to a size of the microstructure e weakly converges in H 1 {D) d as e — > to the 
solution u* of the homogenized problem (4.18). This u* is a unique minimizer of a T— limit 
for the sequence of functionals (-E^) e >o corresponding to the original family of the problems. 
It is also shown in Theorem 4.12 that a first order corrector of u* constructed based on the 
solutions of the cell problems yields a strong H l (D) — convergence. 

To indicate possible future directions of generalization of the considered case we mention 
the study of different types of surface forces with coupled electrostatic or magnetic phenomena 
added. Another important and challenging task is to develop computational techniques to 
analyze the relative importance of electrostatic interactions and surface energy in materials 
with ionized inclusions. Ultimately, such a fundamental understanding of particulate flows 
may be used for the control and optimization of the response in the case of the design of novel 
materials. 

6 Appendix 

Here we prove Theorems 2.4, 2.5 and 2.6. 

6.1 Proof of Theorem 2.4 

Proof. We notice first that according to Remark 2.3, the integral / fi(x,u) ■ (fids makes 

Jd\, 

sense if u belongs to HQ(D) d . Let us assume that {u,p} is a weak solution for the system 
(2.9), so — diver = / weakly in Df. This means that 

/ 2fie(u) : e((fi)dx — I pdiv<p= / / • (fidx, 
Jd } Jd s Jd s 

for every <fi from H^Df) 1 *. Now, if we take a (fi from HQ(D) d , such that e(0) = in D r , we 
obtain 

/ 2fie(u) : e((fi)dx — / pdiv0= / f ■ cf>dx + / an-<pds. 

JD f JD f JDj JdDf 

But 

N 

/ an ■ <fids = / an ■ (fids = > / an ■ (fids, 

JdD f JdD r l=l JdD\, 

by our convention that n represents the unit normal pointing outside the fluid region. In each 
D\. the vector field (fi, having the property that e(<fi) = 0, takes the form 

<fi = M t (x -xi) +m h 

where Mi £ M. is a skew symmetric matrix and mi 6 R is a fixed vector. Thus 

/ an ■ (fids = I an ■ (Mi(x — xi) + mi)ds = 
JdDl JdDL 
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/ y 1 (<rn)i(Mi)ij(x — x{)jds + / an ■ raids = 

JdDl -^ ' " ' JdDl 

/ y2(Mi) ij ((crn)i(x-Xi)j-(an) j (x-Xi)i)ds+ an ■ m^s = 

JdD\. r^ JdD\. 

y^(Mi)jj / (an)i(x-xi)j - (an)j(x - xi)ids + mi- / crnds. 
~Tj JdD\. JdD\, 

Using now the balance of forces (2.5) we have that 

/ ands = / fdx+ / fi(x,u)ds, 

JdD l r JD l r JdD l r 

and using the balance of torques (2.6) we have that 

/ {{an)i(x - xi)j - (trn)j(x - x t )i) ds = 
JdD l r 

/ (fi(x-xi) j -f j (x-xi) i )dx+ / ((fi)i(x,u)(x-xi)j - (fi)j(x,u)(x-xi)i)ds. 

JD' JdDl 



Combining the last three equations we obtain 

d 



/ ancf>ds= Y] fi(Mi)ij(x - xi)jdx + / / • n%idx+ 

JdDl JD ri~j^l JD r 

+ / y2(fi)i(Mi)ij(x- xi)jds+ ft -mids 

JdD\. ?-7± x J dD l r 

= f ■ 4>dx + / f r cj>ds. 

JD l r JdD l r 

Adding over < I < N we get 

A? 



/ an ■ cj)ds = f ■ <pdx + \] / // ' 4>d & 

JdD r JD r f^ JdD l r 



so the variational formulation 



/ 2fie(u) : e(<p)dx — / pdiv(f>= / / • <fidx + / an ■ <pds 

J D f ' JDj JD f JdDf 

becomes, because e((f>) = in D r 

N 

/ 2fie(u) : e{(j>)dx — / pdiv</> = / f-(f>dx + 2~] / fi ■ <pds, 

JD JD f JD J^[ JdD l r 

which is what we wanted to prove. 

Conversely, assuming (2.14) we obtain for all </> from H^Df) 

/ 2fie(u) : e(<j))dx — / pd\v<p= / f ■ <pdx, 
Jd } JD f JD f 
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which gives — div a = f in Df. Using this equation, and taking <fi from Hq(D) with e(<fi) = 
in D r we obtain that 



N N 

V/ fi(x,u)-4>ds = y^ crn-(f>ds, 

U JdD 'r U JdDl r 



which gives us in the same way as in first part the balances of forces and torques. □ 

6.2 Proof of Theorem 2.5 



Proof. We notice first that according to Remark 2.3 the integral / gi(x, v)ds makes sense 

JdD\ 

for v in HQ(D) d . From Theorem 2.4, if {u,p} is a weak solution for the system (2.9), then 



/ 2/xe(w) : e(<f>)dx = \ f ■ 4>dx - Y] / V u g t (x,u) ■ <f>ds, 

Jd Jd JT[ JdD\. 



N 

(6.1) 



for all <j> G V r . So 

N 



E r (u+(f>)-E r (u) = / fie(cf>) : e(c(>)dx+S2 (gi(x,u + <j>) - gi(x,u) + V u gi(x,u) ■ <j>) ds, 

Jd J^[ JdD l r 

which is positive, based on the convexity property of gi (see Remark 2.3). The quadratic 
term / /ue(i>) : e(v)dx makes th functional E r strictly convex, therefore the minimizer is 

JD 
unique. □ 

6.3 Proof of Theorem 2.6 

Proof. Let u be the unique minimizer of E r , and let us show that u satisfies (6.1). For every 
A € R and every <f> £ V r we have that 

E r (u + \<f>) - E r {u) > 0, 

which gives after calculations that 

r f o i \ im\a f -eMj , V^ f 9l{x,u + A0) - gi(x,u) 

hm / 2/j,e(u) : e(cp)dx — I fcpdx + > / — -ds > 0, 

x \oJd Jd i^t^dDl A 

and 

hm / 2„e(u) : e{d>)dx - [ f4>dx + V f 9i(*,u + Xcj>) - gi (x,u) ^ ^ Q 
^SoJd Jd i^lJd Dl r A 

Combining these two we obtain (6.1), which is 

N 

/ 2fxe(u) : e((f))dx = / f-^dx-S^ V u gi(x,u) ■ <f>ds, 

Jd Jd l=1 JdD l r 
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for all 4> £ V r . If is a divergence free vector field from H(j(Df) we obtain 

/ 2fxe{u) : e(cf>)dx = f ■ d>dx, 
Jd Jd 

and from here the existence of p, unique up to an additive constant, such that —diver = / 
weakly in Df (see [33]). As in Theorem 2.4 we obtain 

N N 

JZ / -W u gi(x,u) ■ 4>ds = y~] / an- <j>ds, 

JT[ JdD\. f^ JdD l r 

which gives us the balances of forces and torques. □ 
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